Theoretical method for the study of the excited states of a system 



Ramon Alain Miranda QuintangQ 
Department of Radiochemistry, Higher Institute of Technologies and Applied Sciences, 
Ave. Salvador Allende y Luaces, Quinta de los Molinos, 
Plaza de la Revolucion, POB 6163, 10600, Havana, Cuba. 

A novel, exact, theoretical method for the study of the excited states of a system is presented. It is 
demonstrated how to transform the excited state problem of one Hamiltonian into the ground state 
problem of an auxiliary one. From this, a new exact density functional suitable for excited states is 
constructed. These results make the excited states of a system accessible through all ground state 
theoretical techniques. 

CN \ PACS numbers: 03.65.-w, 71.15.Qe, 31.15.ec, 71. 15. Mb 

i Keywords: Excited states, Variational principle, Density functional theory 

o 

(N ■ 

In spite of the development in the course of the years of many sophisticated and powerful theoretical tools, the 
variational principle of Rayleigh and Ritz (RR) keeps its central role in quantum theory. This principle, besides its 
• many practical applications, has also served as starting point for the development of new theoretical methods, the 
\ density functional theory (DFT) being an example of thisir— . But the classical formulation of the variational principle 
■ is restricted to the study of ground states or to the first excited states having a symmetry which is different from 
the symmetry of the ground state. Given the importance of the analysis of the excited states of a system, a great 
effort has been put in the search for generalizations of the RR principle in ways that can be used to study themir— . 
Here we demonstrate a set of theorems that allow us to apply the original RR principle to excited states. In doing 
so, from a Hamiltonian Hq and with some knowledge of its ground state, a new Hamiltonian H± is built, which has 
the characteristic that its ground state exactly coincides with the first excited state of Ho. A DFT type variational 
principle for the new Hamiltonian is obtained, which enable us to construct a new exact density functional suitable 
^ ■ for the study of excited states. These results make the excited states of a system accessible through all ground state 
O^l techniques. 

Let us start giving some definitions. Let Hq be a Hamiltonian and Eo the set of its eigenvalues. We will denote by 
Vo the set of the eigenvectors of H , if E k G E ,V^ fc will represent the subset of Vn associated with the eigenvalue E k . 
The set {\E k f3)} corresponds to an orthogonal basis of V^ fe where f3 can be interpreted as the eigenvalues of a set of 
observables which, along with the energy, constitute a complete set of observables for the system under study^. Let 
us suppose that Eq — miriEa and that we know a {\Eq/3}}. In this case we can build a Hamiltonian Hi according to: 
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Hi =H a + Kj2\Eol3)(E Q p\ (1) 

@ 

In this expression K is a real constant that will be analyzed later. In full analogy as it has been done before, the sets 
Ei, Vi and Vf fc corresponding to Hi can be defined. 

Now, starting from the definition of Hi, some of its properties will be analyzed. In the following, Halmos'B indicates 
the conclusion of a proof. 

Theorem 1: If E x = min(E a \{E a }) , then: Viji, (ip\Hi\ij) > min{Ei,E + K} 



Proof: Because Hn is an observable we can write^: 



Where £^ |c(£ fc /3)| 2 = 1. 
Then: 



V|^),|V> = Y,< E ^)W)\E k € E ,c(E k (3) = {E k B\i>) (2) 

E k j3 



= <^IW>+^$>(£o/3)| 2 (3) 



Doing \c(E f3)\ 2 = t we will have that: 

VlVMVltfilV) >min T ^{<p\H \cp)+Kt (4) 
The minimum is searched over the set T^, = {\ip) : (Eof3\(p) | 2 = t}. 
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It is then evident that, provided (ip\ip) — 1: minj^ (tp\Ho\tp) = .E^ + E\(l — t) so: 

V|l/>),MW> >e{t) = E 1 + {E Q + K-E 1 )t (5) 

Every linear functional e(t) will reach its lower value in the frontier of its domain of definition, in this case, in t = 
or in t — 1. Checking that, e(0) = E\ and e(l) = Eq + K ends up proving the Theorem. ■ 
Theorem 2: V C Vi 

Proof: G V => = Ef) c(E k P)\E k 0); E k G E ,c(^/3) = (£7*0^) So: 

Fiiv> = ffoi^+^xii^w^) ( 6a ) 

Where (E (3'\E k (3) — SE E k Sj3'/3, being <5jj the Kronecker's delta. 
Let us consider the cases: a,)E k ^ Eq and b)E k = Eq. 

a) E k + Eq^ (E /3'\E k /3) = VjS, /?', then: 

= #o|V>) =^fe|^} IV) G Vi (7) 

Thus, an eigenkct of Hq corresponding to an eigenvalue E k ^ Eq will be as well an eigenket of H\ corresponding 
to the same eigenvalue. 

b) E k =E ^ (E p'\E k p) = {E (3'\E (3) = 5 ff <p, then: 

H x \^) = H Q \iP) +KY,<E p')\Eol3') = Eo^j) + K\i/j) = (E + K)\i/j) => \tf>) G Vi (8) 
P' 

Thus, an eigenket of Hq corresponding to the eigenvalue Eq will be as well an eigenket of H± corresponding to the 
eigenvalue E + K. ■ 

From this proof it follows that (E \{£ }) U {E + K} C E x . 
Theorem 3: K > => (Eo\{£ }) U {E + K} = E 1 
Proof: (by reductio ad absurdum) Let us suppose that: 

3E<£ (E \{£ })U{£o+#} : + \0): H^E) = E\E) (9) 

From Theorem 1 and the RR principle: K > => E ^ 2£ - 
For being 7?o an observable we can write^: 

\E) = Y,c(E k p)\E k p)-E k G Eo, c(E k j3) = {E k f3\E) (10) 

But, according to Theorem 2: VE k (3, E k G Eo, I-E&/3) G Vi then we will have, by hypothesis: 

V£ fc /3,£ fc G E Q ,c(E k f3) = {E k f3\E} = (11) 

The contradiction ends up proving the Theorem. ■ 

The previous results allow us to conclude with the following: 

Corollary: K > => Vi = V , even more, Yf k = Yq" , E k = E k + 6 0k K 

Proof: From the proof of Theorem 2 it follows that V k C V x fc , so it suffices to proof that Y 1 k C V k ■ 

Be \ip) G Vf fc , according to Theorem 3 this implies that E k G (EoU-Eb}) U {E Q + K}. 
For being Hq an observable we can write^: 

|V) = J2 c(EvP)\Ek>P) =Y, C W)W) + J2 c(E k ,f3)\E k ,f3);E k , G Eo,c(E k ,/3) = (E k ,0\1>) (12) 

E k ,p P (E k ,^E k )(3 

In this expression when G Vf we will take fc = 
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From Theorem 2: J2p c(E k j3)\E k /3) G Vf fc , then: 

^ c{E k ,$)\E k ,(3) = 1^-^0(^)1^/3) G vf* (13) 

for being Vf fc a linear subspace. 

Denoting the set of all pure states of the system by W we will have, for being Hi an observabl o 9 ' 10 : 

W = © fc Vf fc (14) 
where @ k stands for the direct sum of the different subspaces. But from (fl"4|^ : 

vf*n£vf fc '={|o>} (is) 

k'^k 

Then: 

|V) =Y,< E kP)W) 6Vf ■ (16) 

P 

Let us analyze the case: 

K>E 1 ~E (17) 

from which we can directly infer that K > 0. 

In virtue of the properties proven before results that, with this selection of K, we have built a Hamiltonian whose 
ground state satisfies E\ = miriEi at the same time that Vf 1 = Wq 1 . That is, the study of the first excited state of 
Hq reduces to the study of the ground state of Hi , for which we can make use of the techniques classically prescribed 
for the study of minimum energy states. 

Among the requirements needed for the construction of Hi , the necessity of having solutions for the ground state 
of Hq that forms an orthogonal basis for this subspace stands out. This, besides been computationally expensive, 
can reminds us the known fact that the RR principle can be applied to the study of excited states under certain 
conditions 11 . In effect, if we work with the expression mmp (ip\Ho\ip), where it has been pointed out that the minimum 
must be searched over the set Po = {|^) : (i(>\<p) — 0V| G V^ }, we will obtain the energy and eingenfunctions 
corresponding to the first excited state of Hq. This procedure is different from the one presented here: in the 
expression of the RR principle corresponding to Hi all N particle functions fulfilling the Pauli principle are included. 
That is, the expression min{i/)\Hi\ip) can be used directly, exactly as the original RR principle mm( - 0|i?o|V')j because 
we are dealing with a ground state. In this new method there is no need to know the set Po. 

The other condition for Hi having the desired properties, namely K > Ei — Eq, is of the outermost importance. 
Even when K > 0, but with K not large enough to fulfill the above condition, the application of the RR principle to 
Hi will end in a ground state with energy Eq + K whose eigenfunctions will be those of the ground state of Hq . In 
order to ensure that (|17l) holds true we must work with K large enough. As a way to check if we have made a right 
choice is enough to determine the min{ip\Hi\ip) , if this value is different from Eq + K the chosen K value is correct. 

Another advantage of being able to apply the RR principle to Hi lies in the possibility of formulate a DFT type 
variational principle for the first excited state of the Hamiltonian if we take, in the spirit of Levy's constrained 
search^: 

Ei = min p E[p] (18) 
Where the minimum is searched over all the A-representable densities, being: 



E[p] = F Ho [p]+ J P(r)v(r)dr (19a) 

FhM = inf^M (r + W + Kj^lEoftiEo^ |V> (19b) 

In the above expressions v(r) stands for the external potential present in Hq, T and W correspond to the kinetic 
energy and particle-particle interaction operators respectively. The subindex in Fjj emphasizes the non-universality 
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of this functional, an unfavorable feature that this method shares with another DFT extensions to excited states^i2~— . 
However, it must be noticed that the non-universal character is respect to the excited state density. The functional 
is in fact an universal bifunctional, depending on the electronic density of the excited and ground state o 15 ! 16 . 

Applications of this scheme within DFT can be expected in two different ways. First, the derivation of formal 
constraints (of the type of coordinate scaling, asymptotic behavior, etc.) for density functionals designed for excited 
states. For this, besides using traditional tools as the virial and Hellmann-Feynman theorems, we can use a novel 
property of this functional. The fact that \/K : K > E\ — Eq E[p] always reach the same minimum value in the same 
density functions will contribute to shed some light in the formal structure of these functionals. Second, practical 
applications in DFT can be derived by means of the optimized potential method^. Here it is not a problem that 
the energy partition for the excited state functionals not needs to be the same that for the ground state functionals, 
since the adiabatic connection can help us to develop new exchange and correlation functionals suitable for excited 
statealsiiil. 

To conclude let us point out that the method presented here can be extended to the study of an arbitrary excited 
state if we have enough information about the stationary states below it. That is, the study of Vjf* , Ek € Eo is 
possible using techniques designed to ground states if {\Ek>/3)}VEk' € Ko,Ek' < Ek are known. 
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